QCD Sum Rules study of meson-baryon sigma terms 
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The pion-baryon sigma terms and the strange-quark condensates of the octet and the decuplet 
baryons are calculated by employing the method of quantum chromodynamics (QCD) sum rules. 
We evaluate the vacuum-to- vacuum transition matrix elements of two baryon interpolating fields in 
an external isoscalar-scalar field and use a Monte Carlo-based approach to systematically analyze 
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I. INTRODUCTION 

The meson-baryon sigma terms are important for 
hadron physics as they provide a measure of chiral- 
symmetry breaking and the scalar quark condensate in- 
side the baryon. In particular, the pion-nucleon and the 
pion-Delta sigma terms have received much attention and 
have been extensively analyzed in many problems (see 
Ref. [ij and references therein) . The sigma terms are re- 
lated to the chiral-symmetry breaking part of the QCD 
Lagrangian 

Cm — —m{uu + dd) — miss — ciui + csusi (1) 

which is expressed in terms of SU(3)-flavor [SU(3)f] sin- 
glet and octet pieces 



Ml = uu + dd + ss, 
Us — uu + dd — 2ss, 



(2) 



with ci = —(2m -I- nig)/?) and Cg — {rUg — to)/3 [m = 
(m,; + mj_)/2 is the average light-quark mass]. The 
strength of the SU(3)i? breaking is controlled by the ma- 
trix elements of the octet piece and can be related to 
the resulting baryon-mass splittings through Gell-Mann- 
Okubo mass formula. The sigma terms, which are defined 
in terms of these matrix elements, can be in turn deduced 
from the STX3)p_pattern and the observed baryon-mass 
differences @, Sljl- 

The sigma terms are equivalent to the values of the 
scalar form factors 

m{B{p')\uu + dd\B{p)) = crB{k)vB{p')vB{p), (3a) 
7fi{B*ip',t')\uu + dd\B*{p,t)) = (3b) 



- <^ iP, t')[g'^''Tis' (fc) + yVi^T(fc)]t^^. {P, t), 
mg{B{p')\ss\B{p)) = aUk)vB{p')vB{p), (3c) 
ms{B*{p',t')\ss\B*{p,t))^ (3d) 

- iP, t'W'^i' [k) + p'^P^Fm^B- {p. t), 

at zero momentum transfer, with B = N , T,, A, and 
B* — A, E*, S*, fl. Here, v{p,t) is the Dirac spinor 
for the spin-1/2 baryon, v^{p,t) is the Rarita-Schwinger 
spin- vector of the spin-3/2 baryon, with the spin pro- 
jection t, k ^ [p' — pY is the momentum transfer, and 
(Tj3(.)(/c), crg(,) (/c) and F^^\k) are the scalar and tensor 
form factors, respectively. The minus sign on the right- 
hand sides (RHS) of (|3bp and (j3dp is conventional like in 
the case of the free Delta Lagrangian. The sigma terms 
are also defined via the Feynman-Hellmann theorem as 



UB' 



TTi-j—^ — rh{B\uu + dd\B), 



q—u,d 



dm„ 



(4a) 



J2 m^^^ = ^m{B*it)\uu + dd\B* {!')), 



drrig 

q—u,a ^ 



m.,^^^ = ms{B\ss\B), 



drris 
dmB" 



-7n,{B*{t)\ss\B*{t')), 



(4b) 
(4c) 

(4d) 



where rng and niB* denote the octet- and the decuplet- 
baryon masses, respectively. 

The matrix elements {B^*^\'ss\B'^*'^) represent the 
strangeness content of the baryons, which can be com- 
bined with the pion-nucleon and the pion-Delta sigma 
terms in order to obtain, e.g., the eta-baryon sigma terms 
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CTr,e(-) = ^{B^*^\m(uu + dd) + 2msSs\B^*^) 



3 

-{oBi') + 2CTg(.)). 



(5) 
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The strange quark condensate of the nucleon is of special 
interest, which is expressed by the ratio 

y^^m^, (6) 

{N\uu + dd\N) 
and related to the pion-nucleon sigma term through 

= aj^°V(l - 2/), (7) 

where cr^^ = 32 MeV is the pion-nucleon sigma term 
obtained from the matrix elements of the octet piece 
in Eq. ([1]) using baryon-mass splittings Hence, a 

discrepancy between cr^'' and the directly observed ctjv 
gives a measure of the strangeness content of the nu- 
cleon. The assumption from the Okubo-Zweig-Iizuka 
(OZI) rule as y = implies a ctat which is significantly 
smaller than expectations based on n-N scattering. The 
resulting puzzle can be solved by considering the possi- 
bility of a non-vanishing ss content in the nucleon. This 
interesting issue has been tackled using various theoreti- 
cal approaches, e.g., the chiral perturbation theory gives 
y ~ 0.21 and lattice QCD gives y ~ 0.36 [1]. 

To our knowledge, there are only a few calculations in 
the literature for the meson- hyperon sigma terms. The 
quark condensates of the baryons have been calculated 
in Ref. [3] by means of a Nambu-Jona-Lasinio (NJL) ap- 
proach to QCD. A chiral model has been used in Ref. [8| 
to evaluate the long-range part of the hyperon scalar form 
factors and the pion-octet-baryon sigma terms. When 
considered in the framework of the octet and the decuplct 
baryons, a determination of the the sigma terms is im- 
portant for understanding the role played by the chiral- 
symmetry breaking in the octet-decuplet mass splittings. 
Moreover, since there is no direct coupling of the pion to 
Lambda baryon, the tt-A sigma term cannot be directly 
determined from experiment. Therefore a theoretical de- 
termination of the TT-A sigma term together with the tt-S 
sigma term is crucial as these terms are related to A-E 
mass splitting. Two of us have recently calculated 01 the 
pion-nucleon and the pion-Delta sigma terms by utilizing 
the external- field QCD sum rules (QCDSR), which are a 
powerful tool to extract qualitative and quantitative in- 
formation about hadron properties [1, [l^, [ll|, ■ In this 
framework, one starts with a correlation function that is 
constructed in terms of hadron interpolating fields. On 
the theoretical side, the correlation function is calculated 
using the Operator Product Expansion (OPE) in the Eu- 
clidian region. This correlation function is matched with 
an Ansatz that is introduced in terms of hadronic degrees 
of freedom on the phenomenological side. The match- 
ing provides a determination of hadronic parameters like 
baryon masses, magnetic moments, coupling constants 
of hadrons, and so on. Our aim in this work is to cal- 
culate the scalar quark condensates of the octet and the 
decuplet baryons and the related sigma terms defined 
in Eq. ([3]), by using the external-field QCD sum rules. 
To determine the value of the sigma terms, we evalu- 
ate the vacuum-to-vacuum transition matrix elements of 



two baryon interpolating fields in an external isoscalar- 
scalar field. For our numerical procedure, we use the 
Monte Carlo-based analysis introduced in Ref. [l^ . This 
method provides a more systematic treatment of uncer- 
tainties in QCDSR. 

Our paper is organized as follows: In Section |TT1 we 
present the formulation of QCDSR and construct the rel- 
evant sum rules. We give the numerical analysis of the 
sum rules in Section Hill Finally, we discuss the results 
and arrive at our conclusions in Section HVl 

II. THE DERIVATION OF THE SUM RULES 

In the external-field QCDSR method, one starts with 
the correlation function of the baryon interpolating fields 
in the presence of an external constant isoscalar-scalar 
field Sq, defined by 

z J d'x e^P-'^ {Q\T[m{x)vBm\Q)s,= 

flB(p) + 5,n|(p) + 0(52), (8a) 
^ I d'^x e^P- (0|r[ry^4*(0)]|0)<,^ = 

[fiB-rip) + Sq [tii,rip) + o{si), (8b) 

where r]B and r^g. are the octet- and the decuplet-baryon 
interpolating fields, which are respectively given as 



m = vn{u s). 




Vn^'nl{u^s), T]^, ^ T]!^,{s ^ d, u ^ s). 

(9b) 

Here a, b, c are the color indices, T denotes transposi- 
tion and C — 17^7°. For the interpolating fields of the 
octet baryons, there are two independent local operators, 
but the ones in Eq. (Pa|) are the optimum choices for the 
lowest-lying positive-parity baryons (see, e.g., Ref [l3| 
for a discussion on negative-parity baryons in QCDSR). 
tleip) and Tls* (p) are the correlation functions when the 
external field is absent and correspond to the functions 
that are used to determine the baryon masses. The sec- 
ond terms in Eqs. ^ represent the linear responses of 
the correlators to a small external scalar field Sq, which 
are computed with an additional term to the QCD La- 
grangian: 

A/: = -S \u[x) u{x) + d{x) d{x)\ - Ss gf [s{x) s(x)] . 

(10) 
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Here, Su = Sd = S {Ss) represent the external scalar 
field and (gf ) is associated with the coupling of the 
external scalar field to the u- and the d- (s-) quark. 
Ilgf,) — ng(.) = Hbi.) (Ilgj.)) denote the correlation 
functions in the existence of the external S (Ss) field. 
The external scalar field contributes to the correlation 
functions in Eq. ([5]) in two ways: first, it directly couples 
to the quark field in the baryon currents and second, it 
modifies the condensates by polarizing the QCD vacuum. 
In the presence of an external scalar field there are no cor- 
relators that break the Lorentz invariance, like (gcr^iyg) 
which appears in the case of an external electromagnetic 
field F'^'^. However, the correlators already existing in 
the vacuum are modified by the external field, viz. 

{qq)s = im) - xS{qq), {ss)s = {ss) - xS{ss), 
{qq)s, = (qq) - x'Ssim), {ss)s, = (ss) - x''Ss{ss), 



{g^qa- ■ Gq) s = (gcqcr ■ Gq) - XcSigc^cr ■ Gq), 
(gcsa- ■ Gs)s = {gcscr ■ Gs) - XGS{gcScr ■ Gs), 
{gcqcT ■ Gq)s, = {gcqer ■ Gq) - xh^s {gcqcr ■ Gq), 
{gcScr ■ Gs)s, = {gcso- ■ Gs) - xhSs{gcS(T ■ Gs), 



(11a) 



(lib) 



where x {= x"" = X'^), X {= X^ = x!^), X" and x" are 
the susceptibilities corresponding to the quark conden- 
sates. Similarly, xg (= Xg = Xg)^ Xg (= Xg = Xg)' 
Xq and Xq denote the susceptibilities corresponding to 
the quark-gluon mixed condensates. Here we explicitly 
assume that the u- and the d- (s-) quark fields couple 
solely to the external field S {Ss}- The quark conden- 
sates get modified in the presence of the external fields 
S and Ss as follows: 



d{qiqi) , d{gcqi(T ■ Gqi) , 

- X {qrqi), — = XG{9cqi(T ■ Gq,), 



dm.j 
for i = j, 

d{qiqi) 
druj 

for i ^ j. 



dnii 



(12a) 



> d{gcq,cr ■ Gq,) -j , - ^ v 
-X-'{qiqi), — ^ XG{9cqiCr ■ Gqi), 



dnii 



(12b) 



where we retain the non-diagonal responses of (qq) {{ss)) 
to external Ss {Sq) field via the susceptibilities x and x", 
and similarly for the quark-gluon mixed condensates via 
Xg and Xq. The coupling of the external scalar field to 
the quark is simply taken as g'^ = gf = 1. 

At the quark level, we have 



T[rjN{x)rj^m O) = 2te'^'^e^''''^'TT{st'\xh.C[ST\x)fCjA^,rsf{x)ri5, 

I Oq V J 

T[v^{x)v^{Q)]\o) ^(o\T[TjN{x)r]^{0)]\o)^ (s^^Ss), 
nVE{x)v^m\o) = (o|r[r;^(a;)77^(0)]|o)^^ (s^^Ss), 



T[^A{x)rj^{0)] 



{Sf {x)^.C[Sf {x)fC^,} j.rsf ix)Yl5 
+ Tr [S2'' {x)j,C[Sf {x)fCj,, } j^rsf {x)r^, - 757^5^' {x)j,C[Sf {x)fC 
xrsf{x)j''j5~l5lf.Sfix)j.C[Sfix)fCj^^S2''{x)ri5) , 



(13a) 
(13b) 
(13c) 

(13d) 



respectively for N, S, S and A, and 



^[<(^)^a(0)] 



2ie' 



'^''^e'^''>'^'(TT{St'\x)^'^C[Sfix)fCr}Sf{x) 



+ 2St''{x)j''C[ST{x)rC^^St'{x)], 



(14a) 
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+ Tr{5r'(a;)7''C[5f (x)]^C7^}5r'(x) + Tr{5r'(x)7'^C[5-'(x)]^C7''}5f (x) 
+ 25-' {x)YC[Sf {x)fC^^sf {x) + 2Sf (x)7^C[5f {x)fC^^S!i^' (x) 
+2Sf' {x)YC[Sf {x)YC^''Sl-' (x)) , 



nvn{x)v^m 



T[v^,{x)r,:^m 



T['nA{x)f]M 



^<5u — > <Ss^ , 



(14b) 

(14c) 
(14d) 



r 



respectively for A, S*, S* and Q. 

To calculate the Wilson coefficients, we need the quark 
propagators in the presence of the external scalar field, 
which are written as 



S,{x)=S'(x)+S'{x). 



(15) 



The first term on the RHS is the part of the propagator 
in the absence of the external field, which is given as 



i [ST ^ (0|T[g«(a;)g-*(0)]|0)o 



1 



27r2x4 

gab 



12 

327r2 

„2 



327r2 2 ^^a;2 



gab^2 



192 



{gcqcr ■ Gq) 



i 6°'''m 



X nioX ln(— X ) + 

48 

X {gcqa--Gq)x + 0{m'^J. 



1 It: 



{qq)x + 



29 X 37r2 
i m„5°'^x'^ 



27 X 32 



(16) 



The second term appears in the existence of the external 
field and is given as 



i[sx^{m'i''{x)q\Q)ms, 



Slab 



29 X 37r2 

gab^q 



327r2 

■x"^ In(-x^) + 



A^fe^cG" a''Mn(-a;2) 



gab 



{qq)x 



gab^2 



27 X 32 



48 

{gcqcr ■ Gq)x 



12 

+ ^^o{9-qcr■Gq)+^^, 



,m. 



^^x''{qq)x 

i mq5''''x'^ ■ 
27 x 32 XGidcqcr ■ Gq)x 



(17) 



when the quark and external field have the same fiavor. 
Otherwise, we have 



mix)tmms,, 



gab^q' gab ^2 

{qq) + ^7^XG{9cq'r-Gq) 



12 



.mq , 



X ■ 



192 



i mqS'^^x'^ 



(18) 



48 ^ 27 X 32 

xXg {9cqo- ■ Gq)x\ + 0{ml , 5^,). 



The analyticity of the correlation function allows us to 
write the phenomenological side of the sum rules in terms 
of a double-dispersion relation of the form 



Jo 



ImnB(p) 



(Sl -p2)(s2 -p2) 



-| poo poo 

Re[n^.r(p) = -^ / / 
^ Jo Jo 



Im[n^.]^'-(p) 

(Sl -p2)(s2 -p2 



dsi ds2, 
(19a) 

■ dsi ds2- 
(19b) 



The ground-state hadron contribution is singled out 
by utilizing the zero-width approximation, where the 
hadronic contributions from the Breit-Wigner form to 
the imaginary part of the correlation function is propor- 
tional to the (5-function: 

ImnB(p) =7r2(5(si - ml)5{s2 - ml) {Q\mW{p)) 

X {B{p)\S{uu + dd)\B{v)) (S(p)|%|0) 

+ 7r25(si - ml)5{s2 - ml) {0\tib\B{p)) 

X {B(p)\S{uu + dd)\£{p)) (f(p)|r?B|0) , 

(20a) 

Imn^(p) =7r2^(si - m%)5{s2 - m|) (0|77b|B(p)) 

X {B{p)\S^-ss\B{jp)) {B{p)\rij,\Q) 

+ 7r25(si - ml)5{s2 - ml) {Q\iib\B{p)) 

X {B{p)\S^-8s\£{p)) (^b)felO). 

(20b) 
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for spin-1/2 baryons. The correlation functions for the 
spin-3/2 baryons are similarly expressed. In the presence 
of external field we have transitions to excited baryon 
states which are denoted by Z. 

The matrix elements of the currents 775 and ?7g. be- 
tween the vacuum and the hadron states are defined as 



(0|7?B|i3(p,t)) = Af5u(p,t), (21a) 
(Olryg.|S*(p,t)) = Az?.w^(p,i), (21b) 

respectively for the octet and the decuplet baryons, 
where Ag and Ag* are the residues. For the spin-3/2 
baryons, we make use of the Rarita-Schwinger spin sum, 
which is 



37-7 



3 mg. 



3 mi.. 



r^''(75-Hmg.), 



{j) + mg. ) 



(22) 



where the slash denotes ^ = "PiiTI^- Inserting Eq. (I^T]) 
into Eq. (PO)) and using the definitions in Eq. ([3]), the 
pole structures of the correlation functions in Eq. (fTQ]) 
are obtained as 



A' 



ArA, 



i> + mg i) + ms 



(23) 



■ ass 



— TTig — 
for spin-1/2 baryons and 



^B' ^6' 5 9. Q^Bf 9pl3 2 _ ^^2 



9 2 
— TOg. 



(24) 



for spin-3/2 baryons. Here, the second terms are asso- 
ciated with the transitions to higher baryon states and 
ag£(.) denote the transition matrix elements. 

We can bring the correlation functions fl^ and 11^, 
into the form 



n^(p) = n|(p2)/ + n|'(p2). 



(25a) 
(25b) 



where the ellipsis represents the Lorentz-Dirac structures 
other than (7^,^ and g^,^]/>. Note that one can obtain 
the sum rules at different Lorentz structures. Here, we 
choose to work with the sum rules at the structures ^ and 
j) fQp t]-^g octet and the decuplet baryons, respectively, 
where the latter are completely contributed by the decu- 
plet baryons with J = \ (see, e.^., Ref. [l3| for details). 



One then expresses the correlation function for the 
octet baryons as a sharp resonance plus a continuum af- 
ter Borel transformation: 

Hg(Af2) ^ (2A|mg-| -I- Cb j 



1 ImHg 

M4 



dso, (26a) 



H^(M^) = 2A^mg^ + C^M^ 



M4 



Im H 



M4 



^ e-^°/^^' dso. (26b) 



Similarly for decuplet baryons, we write 
Hg. (M2) = - 2A^.mg. ^ + Cb' 



M4 



1 /■°° ImHg. _ 



dso, (27a) 



H^.(m2) = - ( 2A|.TOg.^ + C'b, 



m, 
ImH 



M4 



M4 



^e-^°/^^'dso. (27b) 



Here w^{„) and denote the continuum thresholds 

in the existence of the Sq and 5^ fields, respectively. 
Above we have defined H^^,, = Hgj.) = Hg(.,) . We have 
included the single-pole contributions with the factors 
Cbm and Cg(,) , which correspond to transition strengths 
to higher baryon states [second terms on the RHS of ([23]) 
and (|24|) upon Borel transformation]. These transition 
terms are not properly suppressed after the Borel trans- 
formation and should be included on the phenomenolog- 
ical side. 

The QCD sum rules are obtained by matching the 
OPE sides with the hadronic sides and applying the Borel 
transformation. We give the resulting sum rules in Ap- 
pendix [XI where we have defined the quark condensate 
aq — —{2Tr)^ (qq), and the quark-gluon-mixed condensate 
{qgcC ■ Gq) = mQ{qq) with the QCD coupling-constant 
squared — Anas- The flavor-symmetry breaking is 
accounted for by the factor / = {ss)/{qq) and the four- 
quark condensate is parameterized as {{qq)'^) = K{qq)'^. 
The continuum contributions are included via the factors 



E. 



E, 



1+X + ^ 



(28a) 
(28b) 



with X — ^^(.j/M^ and x = (w^(,))^/M^. In the sum 
rules, the third terms on the RHS give the contributions 
that come from the responses of the continuum thresh- 
olds to the external field. Here, S{wg{,))^ and (5(w^(,))^ 
represent the variations of the continuum threshold and 
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the coefHcients are calculated by differentiating the con- 
tinuum parts of the chiral-even octet and decuplet mass 
sum rules with respect to the quark mass. These terms 
are suppressed as compared to the single-pole terms, nev- 
ertheless, should be included on the phenomenological 
side if they are large (see Ref. [l^ for a detailed expla- 
nation of this term). The corrections that come from the 
anomalous dimensions of various operators are included 
with the factors L = log(M^/AQ(^^)/ log(^^/Agp^), 
where /i = 500 MeV is the renormalization scale and 
Aqcd is the QCD scale parameter. A variation of the 
renormalization scale as well as that of the QCD scale 
parameter has little effect on the results. 

III. ANALYSIS OF THE SUM RULES 

We determine the uncertainties in the extracted pa- 
rameters via the Monte Carlo-based analysis introduced 
in Ref. In this analysis, randomly selected, Gaus- 

sianly distributed sets are generated from the uncer- 
tainties in the QCD input parameters. Here we use 
ttq = 0.52 ± 0.05 GeV^, b = (g^G^) = 1.2 ± 0.6 GeV^ 
ml = 0.72 ± 0.08 GeV^, and Aqcd 0.15 ± 0.04 GeV. 
The factorization violation in the four-quark operator is 
searched via the parameter k, where we take k = 2 ± 1 
and 1 < K < 4; here {{qq)"^) > (qq)'^ is assumed via 
the cut-off at 1. The flavor-symmetry breaking param- 
eter and the mass of the strange quark are taken as 
/ = {ss)/{uu) = 0.83 ± 0.05 and m, = 0.11 ± 0.02 GeV, 
respectively (for a discussion on QCD parameters see, 
e.g., Ref. O). 

The value of the susceptibility y can be calculated by 
using the two-point function [l^, EBl 

T(,^)^./.^.e-(0|r[.(.).(.)+.».(.),.(0).(0) 
+ J(0)d(0)]|0), 

(29) 

via the relation 

X{qq) = ^nO). (30) 

The two-point function in Eq. ([29]) at = has been 
studied in chiral perturbation theory [19] with the result 

where /^r = 93 MeV is the pion decay constant and £i and 
£2 are the low-energy constants appearing in the effective 
chiral Lagrangian. A recent analysis of tt-tt scattering 
gives ii = -1.9 ± 0.2_and £2 = 5.25 ± 0.04 Using 
these values of £1 and £2 and taking the quark condensate 
ag = 0.52 ± 0.05 GeV^, we find x = -10 ± 1 GeV"!. 
The susceptibility XG is less certain. It is reasonable to 
assume xg = X: however we adopt a larger uncertainty 



for XG and take XG = — 10 ± 3 GeV ^. As we shall see 
below, the final results are insensitive to a variation in 
XG- 

The susceptibilities x" and Xg related to x 

and xGj respectively, in a strai ght forward way by using 
the three-flavor NJL model [tI. I2H: In NJL model, the 
constituent-quark mass is composed of the current-quark 
mass and a dynamical part (M^) that has a purely non- 
perturbative origin: 

Mi — nii + MP; i = u, d, s, 

(32) 

Mf = -2.g,(g,g,) + M;, 

where rrii is the current-quark mass, gs is the four-quark 
coupling and represents a potentially small contri- 
bution that originates from the six-quark coupling. The 
scalar charge of the constituent quark can be defined by 
using the Feynman-Hellmann theorem as 

Qy = {AU \qjqj\M,) = = h + R^o , (33) 

where we have defined Rij = dMf /drrij with the aid 
of Eq. (jlip. In NJL model, dM[/dmi vanishes in the 
diagonal case and we obtain Ra — ~2 gsdl^^gi) / drrii. 
Using the values of Rij as given in Eq. (4.13) of Ref. 0, 
we find 

Ruu _ Ruu _ X _ 1-14 ^ ^ 1^.. 

R,,~' Rss- rf~ o.A2-'-^' ^^^^ 

which implies that < X- Inserting the central values as 
X = -10 GeV-i and / = 0.83 yields x' - -4-5 GeV"!. 
In our numerical analysis we consider the susceptibility 
values x'* = — 4± 1 GeV^^ by allowing a generous uncer- 
tainty and adopt Xq = X''- 

The non-vanishing values of the non-diagonal suscep- 
tibilities Xj XGj X* and Xg' which the sum rules for 
CT^, CT^ and crji solely depend on, lead to some ss con- 
tent for the nucleon and the Delta, and to uu (dd) con- 
tent for the Omega. Such anomalous quark contents 
are OZI-rule suppressed, therefore we expect that these 
non-diagonal susceptibilities should be very small, if non- 
zero. X and X* can be expressed in terms of a correla- 
tion function as in Eq. (j29p and estimated in chiral per- 
turbation theory. In our analysis, we shall treat these 
susceptibilities as free parameters and adopt the ranges 
-2 < X < -1 GeV-i and -2 < x" < -1 GeV'^ We 
also assume x = Xg and x* = Xg- would like to note 
that these ranges are consistent with those in Ref. [2^ . 
which were taken to reproduce the baryon isospin mass 
splittings. Moreover, since the sum rules for cr^ depend 
solely on the non-diagonal susceptibilities, x and xg are 
actually constrained by the strangeness content of the 
nucleon, which can be determined independently by us- 
ing other approaches. As we shall see below, these values 
of the susceptibilities produce a strangeness content for 
the nucleon in agreement with the expectations based on 
lattice QCD and chiral perturbation theory. 
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TABLE I: The parameter values that we use for the numerical 
analysis of the mass sum rules (|Bl|) - HB8p . the obtained values 
of the overlap amplitudes and the continuum contributions 
for each sum rule at the lower ends of the valid Borel regions. 



Res. 


Region 


cont. 


m (GeV) 


wgiy, (GeV) 


\% (GeV'') 


N 


[0.9-1.2] 


24% 


0.939 


1.5 


1.64 ±0.18 


A 


[0.9-1.4] 


13% 


1.116 


1.7 


2.96 ± 0.30 


E 


[0.8-1.3] 


30% 


1.189 


1.7 


2.90 ± 0.32 




[0.9-1.4] 44% 


1.321 


1.7 


3.76 ± 0.44 


A 


[1.0-1.2] 


30% 


1.47 


1.7 


5.04 ± 1.04 


E* 


[0.9-1.3] 


13% 


1.385 


1.8 


5.36 ± 0.88 




[0.9-1.4] 


7% 


1.533 


2.0 


8.38 ± 1.24 


Q. 


[1.0-1.6] 


6% 


1.672 


2.3 


12.70 ± 1.80 



We use 10'^ such configurations from which the uncer- 
tainty estimates in the extracted parameters are obtained 
using a fit of the LHS of the sum rules to the RHS. For 
TO, we make use of the Gell-Mann-Oakes-Renner rela- 
tion, which is 

2m{qq) = -mlfl (35) 

where = 138 MeV is the pion mass. We use the 
chiral-odd mass sum rules given in Appendix [B] for nor- 
malization of the sigma-term sum rules, which have 
been found to be more reliable than the chiral-even sum 
rules [l^ [2^ . The Monte Carlo analyses of the sum rules 
are performed by first fitting the mass sum rules (jBip - 
(|B8I) to obtain the pole residues Ag and Ag. , and these 
residue values are used in the sigma-term sum rules (|A1[) - 
(|A14[) for each corresponding parameter set. 

The valid Borel regions are determined so that the 
highest-dimensional operator contributes no more than 
about 10% to the OPE side, which gives the lower limit 
on the valid Borel region and ensures OPE convergence. 
The upper limit is determined using a criterion such that 
the continuum-plus-continuum-change and plus-excited- 
statc contributions are less than about 50% of the phe- 
nomenological side, which is imposed so as to warrant 
the pole dominance (this constraint is slightly released 
for a A sum rule |l|). Note that, while the first criterion 
is rather straightforward, one does not initially have a 
complete control on the second, since the phenomcnolog- 
ical parameters are determined from the fit and they are 
correlated. We use the following strategy: we first make 
the fits in a reasonably selected Borel region, which is 
then adjusted by trial and error according to the fit re- 
sults until the above criteria are satisfied. 

The parameter values that we use for the numerical 
analysis of the mass sum rules (|Bll) - (jB8p . together with 
the fitted values of the overlap amplitudes, are given in 
Table [D We also give the continuum contributions for 
each sum rule at the lower ends of the valid Borel re- 
gions (continuum contributions amount to 50% of the to- 
tal phenomenological side at the higher ends of the Borel 
regions) . In order to reduce the uncertainties in the final 



results as much as possible, we fix the baryon masses at 
their experimental values and the continuum thresholds 
at around the first-excited resonance masses as suggested 
by the Particle Data Group [23|- It is well-known that 
the chiral-odd sum rules are less prone to higher order 
corrections in as [25j and they perform better as com- 
pared to the chiral-even sum rules due to cancellations 
in the continuum, however the chiral-odd sum rule some- 
what overestimates the mass of the Delta resonance (see 
Refs. [ll,[2l,[2| for details). Therefore, we take as input 
the value of the Delta mass as suggested by the chiral-odd 
sum rule. 

FIG. 1: The subtracted form of the sum rules ^K^, (|X2|| . 
and (|A3|) . The solid line is the double-pole contribution and 
the dashed-line is the OPE-minus-excited states and minus- 
continuum-change contributions, where we use the average 
values of the QCD and the obtained fit parameters. The error 
bars at the two ends represent the uncertainties in the QCD 
parameters. 
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FIG. 2: (color online) Same as Fig. [T] but for the sum 
rules (|A4|) . (|A5|) and (|A6|) . The error bar at the higher end 
for (jh» is slightly shifted for clear viewing. 
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It is relevant to point out that the dominant contribu- 
tions to the OPE sides of the sum rules (jAip - (jA14[) come 
from the terms that involve the susceptibilities, whereas 
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TABLE 11: The values of the sigma terms, , and the transition strengths, Cgl^-, , for each resonance as obtained from 

a fit of 10^ parameter sets. The second error in o"^',) is due to the uncertainty in the k parameter and the first is the sum 
of errors due to all remaining sources. The second column shows the valid Borel regions and the third column shows the 
continuum-plus-continuum-change and plus-excited-state contributions at the lower ends of the valid Borel regions. We give 
the ratios of the sigma terms in the last column. 



Resonance 


Region 


cont.-|-exc. 


Cb (GeV=) 


QCDSR 


(MeV) 
Rof. [7] 


Ref. [8] 


0"fl/t7A 


A 


[0.9-1.8] 


14% 


1.32(1.15) 


14(03) (02) 


45.4 


33.5 


1 


N 


[0.9-1.3] 


30% 


7.45(2.92) 


53(09)(15) 


56.1 


46.0 


3.73(67) 


E 


[1.1-1.7] 


30% 


7.91(2.96) 


43(07) (08) 


33.3 


29.2 


3.07(62) 


s 


[0.8-1.2] 


28% 


-3.50(1.03) 


7(2) (2) 


27.1 


12.0 


0.44(16) 


Resonance 


Region 


cont.-|-exc. 


Cb' (GeV^) 


(TB' 


(MeV) 






S* 


[1.2-1.6] 


36% 


24.25(6.49) 


56(10)(12) 


23 


- 


1 


A 


[1.3-1.5] 


56% 


39.94(12.65) 


54(10)(15) 


27.6 


- 


0.96(09) 




[1.2-1.7] 


25% 


11.80(2.89) 


28(04) (06) 


17 


- 


0.51(03) 




[1.3-2.1] 


21% 


2.33(1.41) 


6(1)(1) 


9 


- 


0.11(02) 


Resonance 


Region 


cont.-|-exc. 


Ci (GeV""') 




(MeV) 




cbI<^a 


A 


[0.9-1.4] 


21% 


1.95(52) 


243(72)(31) 


191 




1 


N 


[0.9-1.4] 


28% 


0.93(47) 


161(41)(25) 


58 




0.69(19) 


E 


[0.8-1.2] 


30% 


-2.81(39) 


129(36)(34) 


215 




0.58(25) 




[1.0-1.4] 


17% 


-2.44(47) 


272(99) (44) 


341 




1.18(16) 


Resonance 


Region 


cont.-|-exc. 


CIj. (GeV^^) 




(MeV) 






E* 


[1.2-1.6] 


35% 


7.24(2.23) 


301(080)(056) 


138 




1 


A 


[1.2-1.5] 


41% 


4.25(1.66) 


173(053)(034) 


29 




0.57(12) 




[1.2-1.8] 


30% 


10.38(3.19) 


323(091)(048) 


276 




1.09(12) 


n 


[1.1-1.8] 


16% 


11.03(2.91) 


380(117)(048) 


408 




1.29(22) 



the terms that the continuum effects enter with, either do 
not involve the susceptibihty or are proportional to quark 
masses, which are suppressed as compared to the leading 
OPE terms. This leads also to a suppression of the con- 
tinuum contributions. Since the continuum contributions 
in the sum rules (jAip - (jA14|l are suppressed as compared 
to the total phenomenological side, it becomes difficult 
to extract information about the continuum thresholds 
from the fit. Therefore, we have assumed that the con- 
tinuum thresholds are equivalent to those for the mass 
sum rules. We also take w^i,) = The variation 

of the continuum thresholds, S{wg(,))'^ and S{wgf,^)^, 
can also be determined from the fit. However, instead of 
taking these as free parameters, we proceed with a gener- 
ous assumption that the continuum thresholds change by 
25% with the external field viz. (5(u'5(.))^ = Wg{t) /4 and 
)^ ~ /4- We observe that such changes min- 
imally contribute to the final results therefore can safely 
be neglected. 

To demonstrate how well the sum rules and the fitting 
work, we first arrange the sum rules in the subtracted 
form 

ng(.) = A^(.)me(.)— 8(*>' , (36a) 



FIG. 3: Same as Fig.[T]but for the sum rules (|^ and (|A10|I . 
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rris 

where Hgc.) and ng(,) represent the OPE-minus-excited 
state and minus-continuum-change contributions, and 
then we plot the logarithms of both sides. As the RHS 
appears as a straight line with this form, the linearity of 
the LHS gives an indication of the OPE convergence and 
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FIG. 4: Same as Fig. [l]but for the sum rules (|A8|I and 

3 



9j. FIG. 6; Same as Fig.[T]but for the sum rules (|XT2)) and (|XT4)| . 





M"^(GeV"^) 



the quality of the continuum model. This procedure is 
equivalent to searching for a plateau region as a function 
of the Borel mass as in the 'traditional' analysis of the 
QCDSR. Figs. [TJl] show the logarithms of the subtracted 
forms in (|36[) as a function of inverse Borel-mass squared. 
Almost linear behavior of the subtracted forms in these 
figures implies that the valid Borel regions selected ac- 
cording to the criterion above match the plateau regions. 



FIG. 5: Same as Fig.[T]but for the sum rules IfKU) and (KT3\i . 
The error bar at the higher end for is slightly shifted for 
clear viewing. 
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In Table |TT] we present the values of the sigma terms, 
(7g{*) and (7g(t), and the transition strengths, Cg{*) and 
for each resonance as obtained from the fits of 
(|Aip - (IA14|) with lO'^ parameter sets. The second col- 
umn shows the valid Borel regions that are determined 
according to the criterion explained above and the third 
column shows the continuum-plus-continuum-change and 
plus-excited-state contributions at the lower ends of the 
valid Borel regions. For comparison purposes, we also 
give the values of the sigma terms as obtained from a 
chiral model in Ref. 8] and from NJL model 0, IMl (the 
latter are obtained by multiplying the quark condensates 



of the baryons given in Ref. 7] with the central values 
of the quark masses, m = 6 McV and = 110 MeV). 
We quote two errors for the extracted sigma terms: The 
second is the error due to uncertainty in the k parameter 
and the first is the sum of errors due to all remaining 
sources. 

The large errors in the final results and the depen- 
dence on QCD parameters can be substantially removed 
by considering the ratios of the sigma terms. We extract 
these ratios, which are listed in the last column of Ta- 
ble [ni by dividing the corresponding values of the two 
sigma terms for each QCD parameter set and by making 
a statistical analysis of the final distribution. By com- 
paring the values of the sigma terms with those of the 
ratios, we find that the ratios can be determined rather 
accurately due to cancellations in the systematic errors. 
The errors of '^50% in the final results are reduced to 
a level of 10-20% in most of the cases when the ratios 
are considered. We demonstrate this fact by studying 
the correlations between the fit and the QCD parame- 
ters via scatter plots. In Fig. [71 we present the scatter 
plots showing the correlations between ds, tiH*, and x, 
XGi In these figures, the data are normalized with the 
mean values of the sigma terms (the normalized value is 
represented by ag{,) ) so that the results can be compared 
on the same scale. The shaded regions represent the ex- 
tracted values of the sigma terms with their errors. The 
data with the error bars (in blue) are shown for reference 
and give the value of the sigma terms when x, xg or k 
are changed by 3 standard deviations. In Fig. [51 similar 
scatter plots are given for the correlations between a-, 
(t|. , and x*, xh' '^^ "^^^ Monte Carlo analysis has the 
advantage that it covers wide ranges of parameter values. 
We observe no correlation with xo and x^, which implies 
that the results are almost independent of xg a-nd Xg- ^ 
slight positive correlation is observed with the absolute 
value of Xj while the correlations with x" are somewhat 
stronger. The scatter plots showing the correlations be- 
tween the sigma terms and the k parameter imply that 
the sum rules have the strongest dependence on k. This 
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FIG. 7: (color online) Scatter plots showing the correlations between cte, <ys.* , and X) XG, The data are normalized with 
the mean values of the sigma terms (the normalized value is represented by (t^c*) )• The shaded regions represent the extracted 
values of the sigma terms with their errors. The data with the error bars (in blue) are shown for reference and give the value 
of the sigma terms when x, XG or k deviate by 3 standard deviations. 




FIG. 8: (color online) Same as Fig. Qbut for the correlations between a^t , and XGi 




dependence is also the main source of error in the final 
results of the sigma terms. 

We now turn to the ratios of the sigma terms. In Figs. [5] 
and [TOl similar scatter plots are shown for the correlations 
between the ratios a^/a^, a-E^ja^*, cri/at^, a^./ay,' and 
the parameters x, X*, Xg, Xg' impressive to ob- 

serve that the errors are reduced to a great extent when 
the ratios are considered and only a slight dependence 
remains on the QCD parameters. Even the strongest de- 
pendence on K parameter is removed. This behavior is 
common to all sigma terms, which suggests that more ac- 
curate and reliable results can be obtained by considering 
the ratios. 



IV. CONCLUSIONS AND DISCUSSION 

We have derived the QCDSR for the scalar quark con- 
densates of the octet and the decuplet baryons. This, 
together with the definitions in Eqs. ^ and ([5]), leads to 
a determination of meson-baryon sigma terms. We have 
applied a Monte Carlo-based analysis of the sum rules, 
where we have generated randomly-selected parameter 
sets from the uncertainties in the QCD input parameters 
and have made 10'^ fits in order to determine how the 
initial errors propagate to the final fit parameters. To 
determine the valid Borel windows, we have applied two 
criteria which take account of the pole dominance and 
the OPE convergence. 
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FIG. 9: (color online) Same as Fig. [7]but for the correlations between the ratios (Je/o-a, as.* j (ty.' , and Xi XG, 
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FIG. 10: (color online) Same as Fig. [7]but for the correlations between the ratios cr|/crA, o-|* /(Te» , and x' , Xg, k. 
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The large errors in the fit parameters originate from 
several sources. The factorization violation parameter, 
K, and the susceptibilities are the main sources of uncer- 
tainty. We have determined the susceptibihty x of the 
light quark sector in a model independent way by using 
the chiral perturbation theory results. The susceptibil- 
ity can be related to x by using the three-flavor NJL 
model. It is reasonable to assume Xg = X* XG = X 
(with larger uncertainties). Our analysis shows that the 
sigma terms have no considerable dependence on xg and 
Xq: The results are consistent with the current ones even 
for the vanishing values of these susceptibilities. We con- 
sider anomalous quark contents of the octet and the de- 
cuplet baryons by allowing small but non-zero values for 
X: 5C 1 Xg and Xg- Note that such an OZI-rule-violating 
case is considered as a trial complementary to our anal- 
ysis. In the limit of vanishing non-diagonal susceptibili- 



ties, the sigma terms cto, cr^ and cr^ vanish as well (in 
consistency with the OZI rule) and all other sigma terms 
and their ratios are mostly unaffected. 

The overlap amplitudes as determined from the 
baryon-mass sum rules, which we have used to normalize 
the sum rules for the sigma terms, introduce considerable 
uncertainties and can only be improved by a better ac- 
curacy of the mass sum rules as a result of reducing the 
errors in the input QCD parameters. Lastly, the tran- 
sitions to higher-order states and the unknown change 
of the continuum thresholds with the external field are 
another sources of uncertainty. We have observed that 
although a generous range is allowed for the second, it 
has relatively small impact on the final results and can be 
safely neglected. However, the transitions to higher-order 
states, which are not properly suppressed after the Borel 
transformation, should be included on the phenomeno- 
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logical side as they give large contributions for most of 
the cases. 

It is impressive to see that the large systematic er- 
rors cancel when we consider the ratios of the sigma 
terms. We would like to stress the ratio values given 
in Table HIl as our main results, which can be rather ac- 
curately determined and have a minimal dependence on 
the QCD parameters. We have demonstrated this fact 
by studying the correlations between the input parame- 
ters and the sigma terms together with their ratios. We 
have found that the ratios depend very weakly on the 
K parameter as well as on the susceptibilities (within a 
wide range of 3 standard deviations), as a result of which 
the error bars shrink. Our results predict the orderings 

CAT > CTE > CTA > CTH and (TA ~ CTS* > fH* > for the 

pion-baryon sigma terms, independently of the values of 
the susceptibilities and k. As for the strange-quark mass 
contributions to the octet and decuplet baryons, we find 



ctIt > cr^ > cr^ > erf, and (Tq > erf,. ~ (t~, > a 



E ^ii"-! ^ "E* "H* ^ "A- 

The sigma term gives a measure of the contribution of 
explicit chiral-symmetry breaking in the baryon masses. 
The QCD Hamiltonian consists of the chiral-invariant 
terms containing the gauge couplings of gluons and the 
chiral-non-invariant quark-mass term. Suppose that the 
chiral non-invariant term is weak and therefore treated 
perturbatively. Then the sigma term is nothing but the 
contribution of the quark-mass term to the baryon mass. 
We find that, among the octet and decuplet hyperons, the 
chiral-symmetry breaking gives the largest contributions 
to S and S* baryons. A non-trivial outcome is that cth 
and a§2 depend mainly on the non-diagonal responses of 
the quark condensates (besides some small contributions 
from direct couplings). These sigma terms are consistent 
with zero in the limit of vanishing non-diagonal suscepti- 
bilities. We also find that a a is considerably smaller than 
(Te , which indicates that the quark mass term contributes 
to A-E mass splitting. One caveat in the treatment of 
the octet baryons is that instead of the optimum choices 
for the interpolating fields, which are known to perform 
rather successfully in the mass determination 1271 . it is 



possible to adopt a generalized definition in terms of arbi- 
trary mixings between two different local operators. Such 
an extension of our analysis with more general interpo- 
lating fields is desirable, but has the difficulty of treating 
one extra parameter. Note that in the case of the decu- 
plet baryons we have a unique local operator. 

Finally, we would like to comment on the anomalous 
quark content of the baryons. Our sum rules show that a 
non- vanishing response of (qq) to external Ss{ss) field im- 
plies OZI-rule violating ss content for the nucleon and the 
Delta. Similarly, a non-diagonal response of (ss) to exter- 
nal S{uu + dd) field leads to a small but non- negligible uu 
content for the Omega. Therefore, the interesting ques- 
tion of baryon anomalous quark contents boils down to a 
determination of the susceptibilities x, x", xg, and Xg- 
In the range of the values considered for x and Xd the 
value we obtain for the strangeness content of the nu- 
cleon as a% = 161 ± 66 MeV is larger than that from the 
NJL model but compares favorably to expectations based 
on chiral perturbation theory (~ 130 MeV according to 
il), and on lattice QCD (183 ±8 MeV according to 
I). This corresponds to a nucleon strangeness frac- 
tion of y = 0.34 ±0.06 as defined in Eq. M, in agreement 
with that from lattice QCD as y ~ 0.36 Q. The ratio of 
the OZI-rule violating contributions to the nucleon and 
the Delta is of special interest: 




= 0.99 ±0.19, 



(37) 



which predicts an equivalent strangeness content for the 
nucleon and the Delta. 
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APPENDIX A: THE QCD SUM RULES FOR THE MESON-BARYON SIGMA TERMS 

In this Appendix, we give the sum rules for the meson-baryon sigma terms, which are obtained by matching the 
OPE sides with the hadronic sides and applying the Borel transformation: 



CTE : 



^.X al L^/' - ^ a, L-'^/^' + ^{XG + x) ^ L-^/'' 
+ m^x I a, L'^/' + ^XG ml a, f L'^e/^^ 



A|rnE— + Ce + ( ^ - 2m, / a„ ) Swl l'^/^ ^irai:-n,l),M 
m 



(Al) 



I a, if- 2) M^E^-^KX al{2f - 1) L^/^ -^^^fa^M' L^l^ 

+ 'ffa,M' L-W^^ +'^ix + Xg) al{2f - 1) L'^/^^ + !!i^^ + ^^)f L-^^r 

- ^XG ml a, L-2«/2^ - ^xg ml a, f L-^'l^' + ^ (2/ - 1), 
AimA^ + Ca m2 + - ^(3/ - 4) a,) 2.-^9 



y^e'^rnl-wD/M^ 



as : 



Ei - ^a, _ ^ ^ /2 . ^2 ^4/9 



3 



+ ^(X + xg) i-^/^^ + ^XG ml a, f W 
~ m 2 ~ ~ J 



2 



a, (1 + 2/) iJo"^* _ ° ^ (1 + ^2 ^2 ^28/27 _ °^ ^ y ^2 ^2 ^28/27 



- ^(2 + 7/) a, L-l^' + ^(1 + /)(;, + ^2 „2 ^14/27 

+ ^/ (X + xg) i"/^^ + 1|1M6 ijf i'^/^^ - ^^x ^o^* i^/^^ 

+ ^X So^* L'l^-^ + I^XG a, L"^"/^^ 

- ^XG / ml a, L-i"/^^ - 1^ (l + /) « ^1^/27 



xe 



(m=.-«;=.)/M= 



(4/ - 1) Ef L^'l^' -l^nfalM' L^^''' _ ^ /(i + /) a2 ^2 ^28/27 



+ ^(5 - 14/) a, LV27 + 1/ (x + xg) ml a\ L^'l^' 

+ ^/(l + /) (X + xg) ml al L'^'^' + Ef L'^/^^ 

-^xa,M' Ef L'l^' -^fxa, Ef L^/^' + ^Xg a, L-'^'^' 
+ ^XG fa, _ 1^ (1 + /) « 1.16/27 



Al.ms.^ + Ch- + -= 

" m V 5 3 



- 3m, /X a, i;o" L^^^ + ml fa, L-'^l^'^ 

TfX \ O J 



xe 



TTIq ^ 



- 4m, M6 Ef L-^l^ + m,x^ / a, L-^l^ + ri^^s^^laj L-^^'^' 

4m, o 



m. 



^- 2m, /a, 5(^.1)2 M^i-^/s 



xe 



(m=-(«;|)=)/M= 
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^ (4 - 3/) E^-Inxsf a\ L^/' _ 1« (2/ - 1) L^'^ 
+ ^ (/ - 2) a, + ^(X^ + X?;) / al L'^^^ 



TO, 

Is 

4m 



+ ^(X^ + xh) al{2f - 1) L-y-' + 4to, L 



-8/9 



4m, 
9 



m. 



(w^)^ 3m, 
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-m%/M^ 



a£. : 

(4 - /) Ef -^x'^f al £^8/27 _ § ^ + „2 ^28/27 



+ ^(5/ - 14) a, LV27 + + ^2 „2 ^14/ 



27 



m 



4m, 



18 



m 

^g(m|.-(™J.) = )/M 



2m 



2 



^a, (2 + /) Li«/2^ - ^ /(I + /) al ^^8/27 - ^ « / L28/27 

- ^(7 + 2/) a, L^^' + 1/(1 + /) ix' + xh) ml L^'^'' 



9 
7 



2me 



27 



, ^m^ ~s 2 



+ '-^xhmlfa,M'L~''/'' 



m. " 



,s \4 



4m < 



xe 



(m|.-(u)|..)2)/Af=^ 



+ If ix' + xh) ml al L'^'^' - E^ L'^l^' - 3m, Jx' a, E^ L^'^' 

+ ^Xh ml fa, - 4m, f na', L^^''' 



Here M is the Borel mass and we have defined A^(,) — 327r'*Ag(,) and C''^^,^ — l6TT^C^^,y 



APPENDIX B: THE QCD SUM RULES FOR THE BARYON MASSES 



We use the foUowing chiral-odd mass sum rules for normahzation of the sigma-term sum rules [23l |29| . [3 



TUN : 



16 



3 (4 - f)a, - ^(4 - f)a, - ^M^ L-^'^ + ^bM^ E^ L 



-8/9 



+ 



Toe : 



4 TO, 



Kalis - f) 



■ TOA e 



fa,M^Ef - l^+m,M'E^L-^'^ - ^bM^E^L-'/^ + ^na^ 
io o o 



= -TOEe 



-m|/M^ 



Toe 



TO,A : 

-a„E^L''/^'M' - -Eomlagiy^'M^ - -a.bL^^''''' = ^mAe"™-/^', 
Toe* : 

^ (/ + 2) ag * _ 2 + 2)TO^ ag L^^' " ^(/ + 2) i''^'' 



24 



^TOE. e-™--/^' 



TOh» : 

^ (2/ + l)ag £;f * - ^(2/ + l)m2 ag Ef L^^' '^iV + 1) i^'/^' 

+ TO^M^ £;f * - "^M^ Ef + ^1^/27 ^ Ml --Z^', 

TOn : 

^/ ag £f _ 2 _^ ^2 m2 E^ _ ^/ j:i6/2r + E^ 
_ !!^m2 iJo" + 2m, /2 « al L^^/-"' = ^ too e"™"/^' . 



(B2) 



(B3) 

(B4) 

(B5) 



(B6) 



(B7) 



(B8) 
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